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Abstract 

The first observation of the decay B® — ¥ K*°K*° is reported using 35 pb^ 1 of data collected by LHCb in 
I proton-proton collisions at a centre-of-mass energy of 7 TeV. A total of 49. 8± 7.5 B® —> (K + tt~)(K~it + ) 

t— I events are observed within ±50McV/c 2 of the £?° mass and 746MeV/c 2 < m Kn < 1046MeV/c 2 , mostly 

coming from a resonant B Q S — > K*°K*° signal. The branching fraction and the CP-averaged K*° longi- 
tudinal polarization fraction are measured to be B (B° -> K*°K*°) = (2.81 ± 0.46(stat.) ± 0.45(syst.) ± 

/\ 0.34 (fs/fd)) x 10~ 5 and f L = 0.31 ± 0.12(stat.) ± 0.04(syst.). 

PACS: 14.40.Nd, 13.25.Hw, 14.40.Bc 
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1. Introduction 

The decay B° -> K*°K*° is described in the 
Standard Model by loop (penguin) diagrams that 
contain a b — > s transition. The partial width of the 
decay arises from three helicity amplitudes that, 
assuming no additional contributions from physics 
beyond the Standard Model, are determined by the 
chiral structure of the quark operators. Predictions 
obtained within the framework of QCD factoriza- 
tion [1] are B{B° S -> K*°K*°) = (9.1^^ 8 3 ) x 1CT 6 
for the branching fraction and 0.63^q'29 for the 
K*° longitudinal polarization fraction. Predictions 
improve to (7.9±|;|) x 1CT 6 and 0.72^°;^, respec- 
tively, when experimental input is used from B — > 
K*cj) 0[3]. The possibility to use B° -> K* Q K*° 
for precision CP- violation studies to determine the 
phases f3 a and 7 of the CKM matrix [3] has been 
emphasized by several authors [51 |H1 Hj ■ The U- 
spin related channel, B° -» K*°K*°, a b — > d tran- 
sition, has been observed by BaBar [§] , reporting a 
branching fraction of (1.28±g |g =h 0.11) x 1(T 6 and 
h = 0.80^;^ ±0.06 with a signal yield of 33.5t^ 
events. An upper limit for the B Q S K*°K*° 
branching fraction of 1.68 x 10 -3 with 90% con- 
fidence level was reported by the SLD experiment 

ma. 

We present in this Letter the first observation 
of the B® — > K*°K*° decay using pp collisions 
at s/s = 7 TeV at the LHC. The data were col- 
lected during 2010 and corresponds to 35 pb _1 
of integrated luminosity. LHCb has excellent ca- 
pabilities to trigger and reconstruct beauty and 
charm hadrons, and covers the pseudorapidity re- 
gion 2 < 77 < 5. The tracking system consists of a 
21 station, 1-metre long array of silicon strip detec- 
tors placed within 8 mm of the LHC beams. This 
is followed by a four layer silicon strip detector up- 
stream of a 4 Tm dipole magnet. Downstream of 
the magnet are three tracking stations, each com- 
posed of a four-layer silicon strip detector in the 
high occupancy region near the beam pipe, and an 
eight layer straw tube drift chamber composed of 
5 mm diameter straws outside this high occupancy 
region. Overall, the tracking system provides an 
impact parameter (IP) [^resolution of 16 //m + 30 
/im/px( GeV/c), and a momentum resolution <r p /p 
below 8 per mille up to 100 GeV/c. Two ring imag- 



2 The impact parameter is the distance of closest ap- 
proach between a particle's trajectory and its assumed pro- 
duction point. 



ing Cherenkov detectors, one upstream of the mag- 
net, and a second just downstream of the tracking 
stations, together provide a typical kaon identifi- 
cation efficiency of 90%. The pion fake rate, over 
the momentum range from 3 — 100 GeV/c, is be- 
tween 5 and 10 percent. Further downstream is a 
Preshower/Scintillating Pad Detector, an electro- 
magnetic calorimeter, and a hadronic calorimeter. 
The LHCb spectrometer also features a large, five 
station muon system used for triggering on and 
identifying muons. A more detailed description of 
the LHCb detector can be found in |11) . 

To reduce the data rate from the LHC crossing 
rate to about 2 kHz for permanent storage, LHCb 
uses a two- level trigger system. The first level of 
the trigger, implemented in hardware, searches for 
cither a large transverse energy (£r) cluster in the 
calorimeters (£r > 3.6 GeV is a representative 
value during the 2010 run), or a single high trans- 
verse momentum (px) muon or di-muon pair in the 
muon stations. 

Events passing the hardware trigger are read out 
and sent to a large computing farm, where they are 
analyzed using a software-based trigger [12] . The 
first stage of the software trigger relies on the se- 
lection of a single track with IP larger than 125 
pm, pt > 1.8 GeV/c, p > 12.5 GeV/c, along with 
other track quality requirements. Events are subse- 
quently analyzed by a second software stage, where 
the event is searched for 2, 3, or 4-particle vertices 
that are consistent with originating from 6-hadron 
decays. The impact parameter x 2 of the selected 
tracks (IPx 2 ), defined as the difference between the 
X 2 of the primary vertex (PV) built with and with- 
out the considered track, is required to be greater 
than 16 with respect to any PV. The tracks are also 
required to have p > 5 GeV/c and px > 0.5 GeV/c. 
The B® decay vertex must have at least one track 
with px > 1.5 GeV/c, a scalar pt sum of at least 
4 GeV/ c, and a corrected masa^j between 4 and 7 
GeV/c 2 . Additional track and vertex quality cuts 
are also applied. 

Events with large occupancy are slow to re- 
construct and were suppressed by applying global 
event cuts to hadronically triggered decays. These 
included limits on the number of hits in the track- 
ing detectors and scintillating pad detector. 



3 The corrected mass is related to the invariant mass m, 
as m corr = ^/m 2 + \p Tm iss\ 2 + \PTmiss\ , where p Tm iss is 
the missing momentum transverse to the direction. 
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2. Selection procedure and signal yield 

To search for the decay process B® —> 
K* a (K+ir-)K* a (K-ir+) we applied a number of 
offline selection criteria. When a four-track sec- 
ondary vertex is found, the reconstructed momen- 
tum of the B® candidate is used to calculate the 
smallest impact parameter with respect to all pri- 
mary vertices in the event. Tracks are required to 
have pt > 500 MeV/ c, and a large impact parame- 
ter (IPx 2 > 9) with respect to the PV. The differ- 
ence in the natural logarithm of the likelihoods of 
the kaon and pion hypotheses must be greater than 
2 for K + and K~ candidates, and less than for tt + 
and 7T _ candidates. In addition, the K + n~ combi- 
nations^] must form an acceptable quality common 
vertex (x 2 /ndf < 9), where ndf is the number of 
degrees of freedom in the vertex fit) and must have 
an invariant mass within ±150 MeV/c 2 of the nom- 
inal K*° mass (this is around ±3 times its physi- 
cal width @]). The K*° and K*° candidates must 
have pt > 900 MeV/c and the distance of closest 
approach between their trajectories must be less 
than 0.3 mm. The secondary vertex must be well 
fitted (x 2 /ndf< 5). Finally, the _B° candidate mo- 
mentum is required to point to the PV. 

To improve the signal significance, a multivariate 
analysis is used that takes into account the proper- 
ties of the B° s -> K* (K+Tr-)K*°{K-7r+) signal, 
as well as those of the background. It is based on 
a geometrical likelihood (GL) jTHJ, [T3] that uses the 
following set of variables as input: 

• B® candidate impact parameter with respect 
to the closest primary vertex. 

• Decay time of the B® candidate. 

• Minimum impact parameter \ 2 01 the four 
tracks with respect to all primary vertices in 
the event. 

• Distance of closest approach between the two 
K*° trajectories reconstructed from the pion 
and kaon tracks. 

• pt of the B® candidate. 

For a given input sample, the above distributions 
are converted into a set of uncorrelated, Gaussian- 
distributed variables. Two vectors are defined for 



4 This expression refers hereafter to both charge combi- 
nations: K + n~ and K~ 7r+. 




Figure 1: Fit to the K + tt~ K~n + mass distribution of se- 
lected candidates. The fit model (dashed pink curve) in- 
cludes a signal component that has two Gaussian compo- 
nents corresponding to the and B° decays. The back- 
ground is described as an exponential component (dotted 
blue) plus the parametrization indicated in the text (dash- 
dotted green). 

each event indicating its distance to the signal {Si} 
and to the background {Bi} hypotheses by means 
of Xs = S Sf and x% — Yl , where the index i 
runs over the five discriminating variables indicated 
above. The quantity A\ 2 = Xs ~~ x\ ^ s found to 
be a good discriminant between the two hypotheses 
and is used to construct the GL function in such 
a way that it is uniformly distributed in the range 
[0, 1] for signal events and tends to have low values 
for the background. The signal input is generated 
by EvtGen [H] and Pythia 6.4 [16] for the kine- 
matic spectra, and the full detector simulation is 
based on GEANT4 [T7]. 

The GL selection requirement was determined 
by maximizing the signal significance. The GL was 
trained using a fully reconstructed B® — > K*°K*° 
simulation sample for the signal, and a selected 
background sample from the first 2 pb _1 of data, 
which is not used in the analysis. The requirement 
GL>0.24, together with the above selection crite- 
ria, resulted in the mass spectrum in Fig.[l]for the 
selected K + ir~ K~ir + candidates. It is observed 
that the events with masses below the signal re- 
gion have on average slightly higher GL values than 
those with masses above. This indicates the pres- 
ence of a background from partially reconstructed 
B decays. 

To describe the data, we have used a fit in- 
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eluding two Gaussian probability density functions 
(PDFs) centered at the B° and B° masses re- 
spectively, a decreasing exponential and the fol- 
lowing parametrization for partially reconstructed 
B-decays 

AM' ( 1 - — J e(M p -M' )e- fc "' M ®G(M-M'- a } 

(1) 

where is the Heaviside-step function, <£> rep- 
resents the convolution, M' is the variable over 
which the convolution integral is calculated, G(M— 
M'\ dp) is a Gaussian PDF with standard deviation 
<7p and M p and k p are free parameters. The fit re- 
sults are given in Table [l] 

The measured signal yield in a window of 
±50 MeV/c 2 around the B° mass is Nj^+^- k~tt+ — 
49.8 ± 7.5(stat.). The width of the B° peak is in 
good agreement with the LHCb resolution mea- 
sured in decays with similar kinematics such as 
B® —> J/ipcj). The significance of the B® signal was 
determined to be 10.9 a by comparing the log of 
the likelihood between the models with and with- 
out signal. When doing this test, the mass and 
width of the B® and B° mesons were fixed to those 
obtained from independent LHCb measurements of 
B° s -> J/ip<t> and B° -> J/tpK* , respectively. The 
peak at the B° mass, though not significant, is com- 
patible with the B° -> K*°K*° branching fraction 
measured by BaBar [5]. 

As the K*° meson is light compared to the B® 
meson, the invariant masses of the three-body sys- 
tems K + K^tt ± and K + tt^tt ± are rather high, 
above those of the charmed hadrons. This kine- 
matically excludes the possibility of contamination 
from b — > c decays with very short charm flight 
distance, in particular B® — > D~tt + . 

After subtraction of the non-B° component, the 
K + ir~ mass combinations were studied, within a 
±50 MeV/c 2 window of the B® signal, by means of 
a maximum likelihood fit in the (m K +^- , mic-w+) 
plane. Three components are included in the fit, 
namely a double Breit-Wigner distribution describ- 
ing Bg —> x*°K* a production, a symmetrized 
product of a Breit-Wigner and a nonresonant linear 
model adjusted for phase space in the K + tt~ mass, 
and a double nonresonant component. The fit re- 
sult, as shown in Fig.[§ gives (62±18)% K*°K*° 
production. The remainder is the symmetrized 
Breit-Wigner/ nonresonant model. 

The shape of the background mass distribution 
was extracted from a fit to the K + n~ mass spec- 



Table 1: Fitted values of the model parameters for the mass 
spectrum, as described in the text. N s and are the 
number of events for the B® and B° signals, [i s is the fitted 
mass value for the signal and a is the Gaussian width. 
The mass difference between B® and B° was fixed to its 
nominal value [4]. Nf, is the number of background events 
in the full mass range (4900-5800 MeV/c 2 ), and c b is the 
exponential parameter in the fit. M p , a v and k p are the 
'parameters of Eq. |TJ. Finally, f p is the fraction of the 
background associated with Eq. (HI. 



Parameter 


Value 


A^ s 


50.1 ± 7.5 


N d 


11.2 ± 4.3 


H S (MeV/c 2 ) 


5362.5 ± 4.8 


a (MeV/c 2 ) 


21.2 ± 3.3 


N b 


90 ± 10 


c b (lO-^MeV/c 2 )" 1 ) 


-3.37 ± 0.55 


k p (10- 2 (MeV/c 2 )- 1 ) 


5.5 ± 5.3 


f P 


06 + 24 

U - UD - 0.05 


Mp (MeV/c 2 ) 


5170 ± 170 


tip (MeV/c 2 ) 


37 ± 23 



trum observed in two 400 MeV/c 2 wide sidebands 
below and above the B® mass. The number of back- 
ground events to be subtracted was determined 
from the results in Tabled] The sizable K*° con- 
tribution present in this background was taken into 
account. 

A model for B° A"*°F*°(1430), representing 
a broad scalar state interfering with B° — > K*° K*° 
was also studied in the available K + tt~ mass range 
of ±150MeV/c 2 around the K*° mass. The small 
number of events made it impossible to measure 
precisely the size of such a contribution for all val- 
ues of the interfering phase. However, for values of 
the phase away from ir/2 and 37r/2 it was deter- 
mined to be below 12%. Further study of this issue 
requires a larger data sample. 

3. Selection of the control channel 

The branching fraction measurement of B® — > 
K*°K*° is based upon the use of a normalization 
channel with a well measured branching fraction, 
and knowledge of the selection and trigger efficien- 
cies for both the signal and normalization channels. 
We chose B° J/ipK*°, with J/tf> f° r 
this purpose. This decay has a similar topology to 
the signal, allowing the selection cuts to be harmo- 
nized, and it is copiously produced in the LHCb 
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Figure 2: Background subtracted K + w~ and K~ir + com- 
binations for selected candidates within a ±50 MeV/c 2 win- 
dow of the B® mass. The solid blue line shows the projection 
of the 2D fit model described in the text, indicating the K*° 
K*° yield (dashed-dotted red line) and a nonresonant com- 
ponent (blue dotted line), assumed to be a linear function 
times the two-body phase space. The dashed red line indi- 
cates the overall — > K*°X contribution. 



Figure 3: Fit to the mass distribution of selected B° — > 
J/ipK*° events. The dashed red curve is the Gaussian 
component for the B signal. The green dashed-dotted 
line accounts for partially reconstructed B — > J/ipX (see 
Eq. |2j. The pink hatched region accounts for a possible 
B® — > J/ip(j> contamination, parametrized as a sum of two 
Crystal-Ball functions 1201 . The combinatorial background 
is parametrized as an exponential and indicated as a blue 
dotted line. 



acceptance. The presence of two muons in the fi- 
nal state means that B° —> J/ipK*° tends to be 
triggered by a muon rather than a hadron, lead- 
ing to a higher efficiency than for B® — >• K*°K*°. 
The differences in the trigger can be mitigated by 
only considering B° — > J/ipK*° candidates where 
the trigger decision was not allowed to be based on 
muon triggers that use tracks from the decay itself. 

The offline selection criteria for B° — > J/tpK*° 
were designed to mimic those of B® — > K*° K*°. 
In particular, all cuts related to the B® vertex def- 
inition were kept the same. We also used the same 
GL as for the signal. 

The overall detection efficiency was factorized as 
e se 'e* r * 9 . The first factor e sel is the probability of 
the generated tracks being accepted in the LHCb 
angular coverage, reconstructed, and selected. The 
second factor t trig defines the efficiency of the trig- 
ger on the selected events. Both are indicated in 
Table [2j as calculated from Monte Carlo simula- 
tion, along with the number of selected events. 
Note that our measurement depends only on the ra- 
tios of efficiencies between signal and control chan- 
nels. 

The event yield for the selected data was deter- 
mined from a fit to the J/ipK + Tr~ invariant mass 
spectrum as shown in Fig. [3j In this fit, a con- 



strained J/ip mass was used in order to improve 
the B° mass resolution and therefore background 
rejection. A component for the particular back- 
ground source B® — > J/ip4>, with <fi — > K + K~, was 
included in the fit, with a parametrization defined 
from simulation, yielding the result 8±8 events. 
The complete suppression of this background was 
subsequently confirmed using the Armenteros- 
Podolanski [18] plot for the K*° kinematics. The 
fit model also includes a Gaussian signal for the 
B° meson, and a combinatorial background com- 
ponent parameterized with an exponential function 
and an additional component to account for par- 
tially reconstructed B —> J ' / -0 X [19] . This partially 
reconstructed component can be described as 

{_ 1 ^ M — JVf j 2 
c -2<^^> + ^ ifM<,, 

- (2) 
where the parameters /i, k and M are allowed to 
float. The fitted signal according to this model is 
indicated in the third column of Table [2] 

A small fraction of the selected sample contains 
two alternative candidates for the reconstructed 
event, which share three of the particles but dif- 
fer in the fourth one. Those events, which amount 
to 3.8 % (3.7%) in the signal (control) channels, 
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Table 2: Selection and trigger efficiencies obtained from simulation. The observed yield found for the signal and control 
channels in the full mass range are also indicated. The efficiency errors are statistical, derived from the size of the simulated 
samples. 





e sel (%) 


e trig ( % ) 


Yield 


B° -> K*°K*° 
B° -> J/tpK*° 


0.370 ± 0.005 
0.547 ± 0.007 


37.12 ± 0.39 
31.16 ± 0.63 


42.5 ± 6.7 
657 ± 27 


ratio 


0.678 ± 0.013 


1.191 ± 0.027 


0.065 ± 0.011 



were retained for the determination of the branch- 
ing fraction. 

4. Analysis of K*° polarization 

The four-particle K + it~ K~n + angular distribu- 
tion describing the decay of B® into two vector 
mesons (K*° -> K+ir' and K*° -> K-ir+) is de- 
termined by three transversity amplitudes Al, Aii 
and A± . The relative fraction of these can be deter- 
mined from the distribution of the decay products 
in three angles Q\, 82 and p. Here 8\ (62) is the 
K + (K~ ) emission angle with respect to the direc- 
tion opposite to the B® meson momentum in the 
K*° (K*°) rest frame, and ip is the angle between 
the normals to the K*° and K*° decay planes in 
the B® rest frame [5]. We will refer generically to 
the 8 angle from now on, unless differences between 
61 and 6 2 become relevant for the discussion. In a 
time-integrated and flavor-averaged analysis, and 
assuming the B® mixing phase /3 S rs as in the 
Standard Model, the angular distribution is given 
by US El] 



d 3 r 



d cos 6\ d cos 6*2 dip 



\A L \ cos 2 9 1 cos 2 6 2 



-J- |-4ii | 2 \ sin 2 6»i sin 2 8 2 cos 2 ip 
1 L 2 



J— 1-4 J- 1 2 7: sin 2 6»i sin 2 8 2 sin 2 p 
1 H 2 



=- \Al I |^4|| I cos <5|| —^-= sin 20 1 sin 28 2 cos p 

II 2v 2 



We denote the polarization fractions by 



(3) 



fk 



\A L \* + \A\\\* + \Ax\* 



and consequently f L + f\\ + f± = 1. No CP vio- 
lation in the mixing or in the decay has been con- 
sidered. The interference terms related to the A± 
amplitude, both proportional to sin^> s , have been 
neglected. Tl,h are the total widths of the low 
and high mass eigenstates of the B® meson, re- 
spectively, and <5|| is the phase difference between 
Al and An . The total decay width is defined as 
T = (T L +T H )/2 and Ar = T L -T H . Note that as 
a consequence of time integration the relative nor- 
malization acquired by the CP-even and CP-odd 
terms is different. The values Ar = (0.062+°;°^) x 



+ °-ui 2 ) x 10 
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1 m 



, k = L,l±. (4) 



10 12 s" 1 and T = (0.6791^ 
used. 

The detector acceptance is compatible with be- 
ing constant in ip. In contrast, it has a signifi- 
cant dependence on the K*° polarization angle 8. 
The two-dimensional angular acceptance function 
e(cos 8\ ,cos 8 2 ) was studied with a full detector sim- 
ulation. It drops to nearly zero asymmetrically as 
cos 0i,2 becomes close to ±1, as a consequence of 
the minimum p and px of the tracks imposed by 
the reconstruction. 

The Monte Carlo simulation of the K*° accep- 
tance was extensively cross-checked using the B° 
—> J/ip K* Q control channel, taking advantage of 
the fact that the K*° polarization in this chan- 
nel was measured at the B-factory experiments 
[HIES]- The function e(cos 6>i,cos 8 2 ) has been pro- 
jected onto the K*° and K*° axes separately, show- 
ing no appreciable difference, and a small average 
correlation, given the size of the simulated sample. 
We have then used the one-dimensional acceptance 
eg (cos 8) as the basis of our analysis, and deter- 
mined it in five bins of cos 9. Since the longitudinal 
polarization fraction for the B° J/ip K*° channel 
is well measured, a comparison between data and 
simulation is possible. Agreement was found in- 
cluding variations of the angular distribution with 
longitudinal and transverse K*° momentum. In the 
region cos 6* > 0.6 these variations were four times 



larger than for lower values of cos 9. 

The background cos 9 distribution was studied in 
two 200 MeV/c 2 sidebands, defined below and above 
the Bg signal region. Like the signal, it showed a 
dip close to cos 9 = +1 and it was parameterized 
as eg ■ (1 + j3 cos 9). A one parameter fit for /? gives 
the result (3 = -0.18 ± 0.13. 

An unbinned maximum likelihood fit was then 
performed to the data in a ±50 MeV/c 2 window 
around the B® mass, in the region cos 9 < 0.6, ac- 
cording to the PDF 

F(e u 9 2 ,<p) - (1 - a)e e {9 1 )eg{9 2 )I(9 1 ,9 2 ,i P ) 
+a{l + f3 cos + cos e 2 )ee(9i)e e (9 2 ). 

(5) 

The background fraction a was determined from 
the fit to the B® mass spectrum described in Sec. 
[2j Only three parameters were allowed to vary in 
the fit, namely ft,, /m and the phase difference 5\\. 

One-dimensional projections of the fit results are 
shown in Fig. [4] The consistency of the measure- 
ment in various regions of the K*° phase space, and 
of the impact parameter of the daughter particles, 
was checked. The experimental systematic error 
on fi was estimated from the variation of the mea- 
surements amongst those regions to be ±0.03. 

The acceptance for B° s -» K*° K*° is not uni- 
form as a function of proper decay time due to the 
cuts made on the IP of the kaons and pions, and 
a small correction to the polarization fractions, of 
order 3%, was applied in order to take into account 
this effect. It was calculated from the variation in 
the measured polarization amplitudes induced by 
including a parametrization of the time acceptance 
in Eq. [5] Note the different correction sign for each 
polarization fraction, as a consequence of the as- 
sumption Ar 7^ 0. 

The sensitivity of the Jl measurement with re- 
spect to small variations of the cos 9 distribution 
has been tested. These variations could be at- 
tributed to experimental errors not accounted for in 
the simulation or to interference with other partial 
waves in the Kn system. A high statistics study us- 
ing B° — > J/tp K*° muon triggers revealed a small 
systematic difference between data and simulation 
in eg (cos 9) as cos 9 approaches +1, which was taken 
into account as a correction in our analysis. When 
this correction in varied by ±100%, fi, varies by 
±0.02 which we consider as an additional source 
of systematic error. The total systematic on Jl is 
thus ±0.04. 
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Figure 4: cos 9 (above) and ip (below) acceptance corrected 
distributions for events in the narrow window around the 
B® mass. The blue line is the projection of the fit model 
given by Eq. [3] for the measured values of the parameters 
/i, /|| and fin. The dotted lines indicate ±1<t variation of 
the fi, central value. 

We finally measure the K*° longitudinal polar- 
ization fraction f L = 0.31±0.12(stat.)±0.04(syst.), 
as well as the transverse components /m and f±. In 
the small sample available, the CP-odd component 
f± appears to be sizable f± = 0.38 ± 0.11(stat.) ± 
0.04(syst.). A significant measurement of Sn could 
not be achieved (Su = 1.47 ± 1.85). 

As seen in Eq. d3j) , due to a nonzero Ar time in- 
tegration changes the relative proportion between 
the various terms of the angular distribution, with 
respect to their values at t = 0. If we call the 
polarization fractions we would have measured un- 
der the assumption Ar = 0, it can be derived from 
Eq. [3] that our measured values are 

fk = fi(l+Vk^f) (6) 

with CP eigenvalue rjk = +1, +1, — 1 for k = L,\\ 
,±. Given the current knowledge of AT/T [4], the 
magnitude of the correction to fk amounts to 4.6%, 
and the associated systematic error related to Ar 
error is 2.6%, which we have neglected in compar- 
ison to other sources. 

5. Determination of the branching fraction 

The results of the previous sections can be 
brought together to provide a determination of the 
branching fraction of the B Q S -» K*°K*° decay 
based upon the use of the normalization channel 
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B° -> J/ipK*° through the expression 



B (B° s -> K*°K*°) = X fL 



trig 

X — 7 — : X 



c sel 

B°^K»°K»° 
N B«^K*°K*° 



trig 
e BO^,K*°K*o 



, N B n^j/Tp K ,a 



xB vis (B° J/#T°) x ^ x H (7) 

where B V i S (B° —> J/ipK* ), the visible branching 
ratio, is the product 6(5° -> J/ijK* ) x B(J/0 -> 
H + x £>(if*° — > if + 7r~). The numerical value 
of 23(B° -> J/tJjK*°) = (1.33 ±0.06) x 10~ 3 is taken 
from the world average in [5J, B(J/ip — > ,u + /i~) = 
0.0593± 0.0006 4] and S(if*° -> #+7T-) =2/3 0. 
The ratio of &-quark hadronization factors that ac- 
counts for the different production rate of 5° and 
B° s mesons is f s /f d = 0.253 ± 0.031 [24 . The fac- 
tor 9/4 is the inverse square of the 2/3 branch- 
ing fraction of K*° —> K + tt~. The number of 
candidate events in the signal and control channel 
data samples are designated by N BO _y K , g*o and 

Nbv^J/^K* - 

The correction factor A f L is motivated by the fact 
that the overall efficiency of the LHCb detector is a 
linear function of the K*° longitudinal polarization 
fi, . Taking into account the measured value and er- 
rors reported in section |4j Monte Carlo simulation 
was used to estimate Xf L = 0.812 ± 0.059. 

We have considered two sources of systematic 
uncertainty associated to the ratio of selection ef- 
ficiencies. The first source results from discrepan- 
cies between data and simulation in the variables 
related to track and vertex quality, and the sec- 
ond is related to particle identification. A small 
difference observed in the average impact parame- 
ter of the particles was corrected for by introducing 
an additional smearing to the track parameters in 
the simulation [23]. While the absolute efficien- 
cies vary significantly as a function of vertex reso- 
lution, the ratio of efficiencies remains stable. We 
have assigned a 2% uncertainty to the ratio, af- 
ter comparison between simulation and the B° — > 
J/ip K*° data. The K/ir identification efficiency 
was determined using a sample of B° — > J/ip K*° 
events selected without making use of the RICH 
detectors. As the signal channel contains one more 
kaon than the control channel, a correction factor of 
1.098±0.019 was applied to the branching fraction, 
and a 2% error was assigned to it. The efficiency of 
muon identification agrees with simulation within 



1.1% 26 . All these factors are combined to pro- 
duce an overall systematic uncertainty of 3.4% in 
the ratio of selection efficiencies. The uncertainty 
in the background model in the B® mass fit (±2 
events) contributes an additional systematic error 
of 4.7%. 

Trigger efficiencies can be determined, for partic- 
ular trigger paths in LHCb, using the data driven 
algorithm described in [26] , This algorithm could 
be applied for the specific hadronic triggers used for 
B° -> J/ip K*°, but not for the small B° -> K*° 
K*° signal. The efficiency related to cuts on global 
event properties, applied during the 2010 data tak- 
ing, is determined from J ftp minimum bias trig- 
gers [2S] . The result indicates a trigger efficiency of 
(26.8 ±3.8)%, smaller than the simulation result of 
(31.16 ± 0.63)% shown in Table [2] Although these 
are consistent within uncertainties, we nonetheless 
apply a —9% correction to the ratio of trigger effi- 
ciencies between B a -t J/ijj K*° and B° -> K*°K*° 
channels, taking into account correlations in the 
trigger probability. A systematic error of 11% was 
assigned to uncertainty on the trigger efficiency, en- 
tirely limited by statistics, both in the signal and 
control channels. Detector occupancies, estimated 
by the average number of reconstructed tracks, are 
larger by 10% in the data than in the simulation. 
This implies an additional correction of +4.5% to 
the ratio of efficiencies, since the control channel is 
observed to be more sensitive to occupancy than 
the signal channel. 

An ~ 8% S-wave contribution under the K*° res- 
onance in the B a —> J/ip K*° channel has been 
observed by BaBar [23], and the data in a ±70 
McV/c 2 mass interval around the K*° mass [27] 
yields a (9.0 ± 3.6)% extrapolation to the ±150 
MeV/c mass window. The S-wave background 
doubles for the K*° K*° final state, and it may 
certainly have a different coupling for both chan- 
nels. Our direct measurement reported in section [2] 
of (19±9)% is still lacking precision to be used for 
this purpose. When evaluating the branching frac- 
tion, we have assumed a 9% S-wave contribution, 
and assigned a systematic error of 50% to this hy- 
pothesis. A summary of the various contributions 
to the systematic error can be seen in Table [3j 

Our final result is 

B(B° S -> K*°K*°) = (2.81 ± 0.46 (stat.) 

± 0.45 (syst.) 
±0.34 (/.//„)) x 10" 5 . 
As we have seen at the end of section [4] un- 



Table 3: Estimated systematic error sources in the 
B (Bj — > K*°K*°) measurement. 



Systematic effect 


Error (%) 


Trigger efficiency 


11.0 


Global angular acceptance 


7.2 


S-wave fraction 


5.0 


Background subtraction 


4.7 


P° J/tp K*° and J/tp -> fjtfx 


4.6 


BR uncertainty 


Selection efficiency 


3.4 


Total 


15.9 



equal normalization factors arise upon time inte- 
gration of individual polarization amplitudes with 
well-defined CP-eigenvalues. This has the inter- 
esting implication that the time-integrated flavor- 
averaged branching fraction (Pi) as determined 
above cannot be directly compared with theoret- 
ical predictions solely formulated in terms of the 
decay amplitudes Aj} + A\\ 2 + Aj_ 2 (B ). Meson 
oscillation needs to be taken into account, since 
two distinct particles with different lifetimes are 
involved. Owing to the fact that A± is CP-odd, 
the relationship between these quantities reads as 
follows 

B =B 1 (l + ^(f L +f ll -f x )Y (8) 

According to our measurements of /z, + /|| — fx> the 
correction is small (3% if current values are taken 
for Ar), and we do not apply it to our measure- 
ment. 

6. Conclusion 

The b ->• s penguin decay B a s -> K*°K*° 
has been observed for the first time. Using 35 
pb^ 1 of pp collisions at 7 TeV centre-of-mass en- 
ergy, LHCb has found 49.8 ± 7.5 signal events 
in the mass interval ±50MeV/c 2 around the B® 
mass. Analysis of the K + tt~ mass distributions 
shows that most of the signal comes from B a s — > 
K*° K*°, with some S-wave contribution. The 
branching fraction has been measured, with the 
result B (B° ->■ K*°K*°) = (2.81 ± 0.46(stat.) ± 
0.45(syst.)±0.34 (f s /f d )) x 10~ 5 . The CP-averaged 
longitudinal A'* polarization fraction has also been 
measured to be f L = 0.31±0.12(stat.)±0.04(syst.), 
as well as the CP-odd component f± = 0.38 ± 
0.11(stat.) ±0.04(syst.). 



When we consider our measurement in associa- 
tion with that of [S], it is remarkable that the lon- 
gitudinal polarization of the K*° mesons seems to 
be quite different between P° -> K*° K*° {f L _= 
0.31±0.12(stat.)±0.04(syst.)) and P° -> K*° K*° 
{f L = 0.80tg;^(stat.) ± 0.06(syst.)), despite the 
fact that the two decays are related by a U-spin ro- 
tation. However, the ratio of the branching ratios 
of B° s and B° decays is consistent with 1/A 2 where 
A is the Wolfenstein parameter, as expected. 
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